A switchable controlled-NOT gate in a spin-chain NMR quantum computer 
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A method of switching a controlled-NOT gate in a solid-stae NMR quantum computer is pre- 
sented. Qubits of I = 1/2 nuclear spins are placed periodically along a quantum spin chain (1-D 
antiferromagnet) having a singlet ground state with a finite spin gap to the lowest excited state 
caused by some quantum effect. Irradiation of a microwave tuned to the spin gap energy excites 
a packet of triplet magnons at a specific part of the chain where control and target qubits are in- 
volved. The packet switches on the Suhl-Nakamura interaction between the qubits, which serves 
as a controlled NOT gate. The qubit initialization is achieved by a qubit initializer consisting of 
semiconducting sheets attached to the spin chain, where spin polarizations created by the optical 
pumping method in the semiconductors are transferred to the spin chain. The scheme allows us to 
separate the initialization process from the computation, so that one can optimize the computation 
part without being restricted by the initialization scheme, which provides us with a wide selection 
of materials for a quantum computer. 

PACS numbers: 76.60.-k, 03.67.Lx, 75.45.+J 
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I. INTRODUCTION 

A quantum computer (QC) is a Turing machine which 
performs information processing based on the principles 
of quantum mechanics. It makes the best use of the fea- 
tures of the quantum mechanics, such as superposition 
and entanglement of quantum states. These features en- 
able us to perform parallel computation for all the pos- 
sible states simultaneously, which makes it possible to 
deal with the problems that are formidable for classical 
(conventional) computers. New algorithms recently dis- 
covered jj], ||, D have shown the great promise of the 
QC's, which has accelerated the attempts to implement 
the QC's in actual physical systems. 

A QC is composed of a set of two-level systems called 
qubits. The qubits need to be isolated enough from the 
environment, and controllable from outside to deal with 
the information. In this respect, nuclear spin systems in 
a matter are the promising candidates, because they are 
only weekly coupled with the environment (electron sys- 
tems) through hyperfine couplings, and controllable by 
the well-established technique of the nuclear magnetic 
resonance (NMR). Actually, the first 2-qubit QC's were 
implemented by solution NMR j|, ||. The quite suc- 
cessful implementations proved the great promise of the 
NMR-QC's. It is unfortunate, however, that the solution 
NMR-QC has a difficulty in its scalability. The number 
of available qubits in solutions is limited because of the 
limited number of nuclei in one molecule. The required 
number of qubits with which QC can surpass its classical 
analogue is estimated to be more than 10 3 , so that it is 
a primary concern to increase the number of qubits. 

So far, two models have been proposed to systemat- 
ically increase the number of qubits in the NMR-QC's. 
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Kane proposed a multi-qubit NMR QC model, which uti- 
lizes 31 P nuclei embedded in a Si matrix 0. This model 
provides us with a chance to extend the number of qubits 
systematically, whereas it requires fine structures to be 
fabricated and controlled, realization of which is a chal- 
lenge to the current state-of-the-art nano-technology. On 
the other hand, Yamaguchi and Yamamoto proposed to 
utilize nuclei as they arc in nature, i.e., in crystals [Q. 
The 31 P nuclei in CeP are placed in the magnetic field 
gradient, so that each nucleus (qubit) can be accessed by 
adjusting the NMR frequency. The proposal is very at- 
tractive because of its simple structure, although it still 
has some technical obstacles to be cleared ||. We have 
pursued this possibility checking the key issues. 

One of the key issues is how to provide an inter-qubit 
(inter-nuclear) coupling, which is used for a controlled 
NOT (c-NOT) gate shown in Fig. [I]. So far, a nuclear 
dipole (direct) coupling has been supposed to be a po- 
tential inter-qubit coupling g|. It is unfortunate, how- 
ever, that the dipolar coupling is always present whenever 
qubits are put close to each other, so that one should con- 
tinue applying the decoupling sequences to remove un- 
wanted couplings. Apparently, these sequences consume 
a great deal of time, causing a longer computation time. 
Moreover, as the number of qubits is increased, the num- 
ber of inter-nuclear couplings is increased accordingly, 
so that the decouplings become more and more com- 
plicated, and they eventually become formidable. So, 
a decoupling-free QC is highly desired, i.e., the inter- 
nuclear coupling should be able to be switched on when, 
and only when necessary. 

Another key issue is the method of distinguishing, in 
the frequency domain, each nucleus (qubit) separated by 
a distance of the order of a lattice constant. Even the 
field gradient of lT//im JTo[ is still marginal for this pur- 
pose |9(]. The larger inter-nuclcar (qubit) distance in the 
real space makes the distance between the adjacent NMR 
lines wider in the frequency domain, which relieves this 
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FIG. 1: Schematic view of the sequence for the c-NOT gate 
in the rotating frame of the target qubit. Upper and lower 
figures correspond to the cases of the up (0) and the down 
(1) states of the control qubit. (a) The target qubit spin is 
pointing to the Z-direction. The spin is rotated by a n/2 pulse 
around the negative X axis, (b) The spin starts to turn in the 
XY plane due to the additional field caused by the control 
qubit. (c) The spin turns by ±90° to positive or negative X 
direction in the XY plane, (d) The spin is rotated again by 
a 7r/2 pulse around the positive Y axis. The direction of the 
spin with respect to Z is hereby controlled according to the 
spin state of the control qubit. 



constraint. Hence, the inter-nuclear coupling should be 
able to reach rather long distance. Unfortunately, the di- 
rect nuclear dipole coupling reaches at most a few lattice 
points jy}. 

These facts motivate us to seek for the possibility of the 
long range indirect couplings mediated by electrons. The 
indirect couplings include the J-coupling due to the cova- 
lent bondings, the RKKY interaction in metals |12[|, and 
the Suhl-Nakamura (SN) interaction in magnets Jl3| , [l4j . 
Among them, the SN interaction has the characteristics 
preferable for the present purpose, such as the long-range 
nature of the coupling and the external controllability of 
the coupling strength. 

In this paper, we present details of the model of a solid- 
state NMR-QC with an inter-qubit coupling provided by 
the SN interaction [Tq] . A singlet-triplet transition in 
a quantum spin chain (ID antiferromagnet) provides us 
with a switch for the inter-qubit coupling. The model has 
the following advantages over the other existing models. 

(1) It is intuitive because the computation starts in the 
silent environments rather than the turbulence of interac- 
tions, which makes the designs of the logic gates simpler. 

(2) The long-range nature of the SN interaction allows 
us to place qubits apart farther from one another, which 
facilitates distinguishment of the qubits in the frequency 
domain. (3) The presence of the spin gap in the quantum 
spin chains allows simple operation of the gate switching 
compared to that proposed for the Kane's type NMR-QC 
p6[ . We also present the scheme of qubit initializer, an 
effective nuclear polarizer comprising the optical pump- 
ing and the polarization transfer methods. The scheme 
enables us to separate the materials responsible for the 
initialization and computation, so that one can optimize 
the computation part without being restricted by the ini- 
tialization scheme. 
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FIG. 2: Schematic illustration of the c-NOT-gate switching 
in a ID antiferromagnet in a magnetic field gradient. The 
1 = and 1/2 (= qubits) nuclei are shown by balls and arrows, 
respectively. (1) All the inter-nuclear couplings are switched 
off in the absence of magnons. A 7r/2 pulse is applied to the 
target (In+i) qubit. (2) A magnon packet (hatched part) 
is excited between the control (In) and the target (7jv+i) 
qubits. The magnon packet entangles the two qubits, and the 
second 7r/2 pulse makes the 7jv+i spin turn either back to the 
+1/2 state for In = 1/2 (2a), or forth to the -1/2 state for 
In = -1/2 (2b) (see Fig. §). The (2a) and (2b) states are 
superposed in the actual computation. 



II. MODEL 

Here, we present the outline of the model. The main 
idea is illustrated in Fig. ||. The system consists of a one- 
dimensional array of electron spins (quantum spin chain) 
placed in a magnetic field gradient, which is produced 
by a micromagnet fabricated outside the spin chain ||To| . 
Suppose that the electron spins are paired into singlets 
in the ground state (| ss z ) = \ 00)) with a finite gap to the 
lowest triplet branch ( | 1 — 1) ) of the k = magnon modes 
because of some quantum effects. Examples of such sit- 
uations can be found in spin ladder, Haldanc, dimer and 
spin-Peierls systems. Also suppose that nuclei (1 = 1/2) 
serving as qubits can be placed periodically, e.g., every 
ten lattice points, each of which has a hyperfme coupling 
with the electron spins. 

Since there are no unpaired electron spins in the 
ground state, the nuclear spins are well-isolated from 
the environment. The nuclei are also decoupled from 
the charge and lattice properties of the electrons because 
/ = 1/2. The rather long distance between qubit nu- 
clei is effective both to diminish the direct nuclear dipole 
couplings between qubits and to distinguish one qubit 
from another in the NMR frequency domain, because 
the longer inter-qubit distance leads to larger interval 
between NMR lines in the given field gradient. 

It is known that virtual exchanges of magnons between 
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nuclear spins result in the inter-qubit interaction called 
the Suhl-Nakamura (SN) interaction Jl3], 14 . Since it 
is not necessarily accompanied by actual excitations of 
magnons, it can exist even at temperatures well beneath 
the excitation gap. Nevertheless, the transverse compo- 
nent of the SN interaction with a form of if I J is always 
absent in the field gradient because of the detuning ef- 
fect, i.e., a mismatch in the Zeeman energies prevents the 
nuclei from exchanging magnons. On the other hand, the 
longitudinal component (ZfiJ) of the SN interaction can 
survive even at low temperatures, which is characterized 
by the antiferromagnctic spin-spin correlation function in 
the ground state (see 
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however, is short-ranged because the spin-spin correla- 
tion function decays exponentially as a function of the 
distance, so that the interaction can reach only a few 
lattice points [^0|. As a result, the interactions between 
nuclei apart from each other by the order of ten lattice 
points are completely switched-off at low temperatures. 

In order for the system to work as a QC, one has to 
provide logic gates. A QC is complete if it is equipped 
with arbitrary rotation (R) and controlled-NOT (c-NOT) 
gates. ||l| The R gate is a single qubit operation, which 
can be accomplished by an NMR pulse with an appropri- 
ate pulse width at the corresponding frequency. On the 
other hand, the c-NOT gate is a two-qubit operation, 
which works in such a way that a target qubit changes 
its logic according to the state of the control qubit, which 
can be realized by the pulse sequence shown in Fig. [I]. 
Since the c-NOT gate between any combination of qubits 
apart from each other along the chain can be realized by 
a series of c-NOT gates between adjacent qubits |22| , the 
QC is complete with the inter-nuclear coupling with the 
form Ii Z Ij Z between adjacent nuclear spins. This cou- 
pling is provided by the longitudinal component of the 
SN interaction via the k = magnon mode composed of 
the spin triplet state p3| , which is selectively excited at 
a specific part of the chain where the c-NOT gate is to 
be performed. 

The transition between singlet and triplet states by a 
microwave irradiation can be used to create the triplet 
k = magnons. Although the excitation is primarily 
forbidden for the usual electric dipolar transition, it of- 
ten becomes possible in the actual systems because of 
some higher order terms in the electron-photon interac- 
tion Hamiltonians [^4], ^5|. The position of the excited 
magnons along the chain can be specified by the applied 
microwave frequency, which is uniquely given in the field 
gradient. The energy diagrams of the k = magnon ex- 
citations in the magnetic field is shown in Fig. || ^1], ^(| . 
In the field gradient, the magnetic field at each part of 
the chain is unique so that the excitation energy to the 
lowest triplet state (| 1 — 1)) is also uniquely given, which 
provides us with a spatial resolution of the excitation, i.e., 
one can specify the position of the magnon excitations by 
adjusting the frequency of the applied microwave. The 
region and the amplitude of the excited magnons can be 
controlled by the frequency resolution and the power of 




IT (100 qubits) Magnetic Field (T) 



FIG. 3: Energy diagrams for the three triplet branches (| 11), | 
10) and | 1 — 1)) of the k — magnon excitation modes against 
the singlet state (| 00)) as a function of the external field. In 
the field gradient, the horizontal axis also corresponds to the 
position along the chain. The thick line on the | 1 — 1) branch 
shows a part of the chain used as a QC, and an open circle is 
the spot where the transition (| 00) — >| 1 — 1)) occurs. 



the microwave. 

Since the two upper states (| 10), | 11)) irrelevant to 
the transition of interest can be ignored, one can intro- 
duce effective spins within the | 00) and | 01) subspace 
p6| . These spins are rotated by the microwave, creating 
a packet of superpositions of | 00) and | 1 — 1) along the 
chain [ p4[ , corresponding to the magnon excitations with 
the wave number k ~ 0. The packet is localized along the 
chain due to the magnetic field gradient (see §111). The 
number of excited spins is given by the balance between 
the excitation to the | 1 — 1) state and the relaxation 
(with the lifetime of T s ) to the ground state. 

The triplet states make an additional field (H tI ) asso- 
ciated with the "shift" at the qubit site, which should be 
distinguished from the additional field from the adjacent 
qubit via the SN interaction (Hsn)- Assuming the hy- 
perfine coupling An = 100 kOe//is and the rotation of 
the spins: n(0) ~ 1 % , H tI is estimated to be ~ 1 kOe, 
which corresponds to ~ 1 % shift of the NMR frequency 
in the magnetic field of 10 T. In practice, the exact value 
of Htr can be measured by the following method; under 
the microwave irradiation, one observes the shift of the 
NMR frequency of the target qubit while saturating the 
control qubit by applying the corresponding NMR rf field 
continuously. The saturation of the control qubit results 
in i?sN = at the target qubit, so that the observed shift 
directly corresponds to H tr . 

Then, the c-NOT gate is achieved as follows (see Fig. 
|l|). (a) In the beginning, the target qubit (Tjv+i) is as- 
sumed to be pointing to the Z direction. The 7r/2 pulse 
with the rf frequency to — j n H(xN+i) is applied in the 
negative X direction in the rotating frame of the target 
qubit. (b) A microwave is applied to the chain. The SN 
interaction is switched on, and the target qubit starts 
to rotate in the XY plane of the rotating frame with 
lj = "f n (H(xis!+i) + H tr )- The direction of the rotation 
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FIG. 4: A spin ladder configuration. Closed and open circles 
represent electron and nuclear spins, respectively. 



with respect to the frame depends on the state of the 
control qubit (In)- (c) By the time t = ir/(2-f n HsN), the 
target qubit reaches either positive or negative X. (d) The 
SN interaction is shut off by shutting off the microwave, 
and the second n/2 pulse with j n H(xN+i) is applied in 
the negative Y direction, which rotates the target qubit 
upward or downward according to the control qubit. This 
sequence allows us to perform the c-NOT gate. 



III. SINGLET-TRIPLET EXCITATIONS IN THE 
QUANTUM SPIN CHAIN 

One of the key phenomena in the present model is the 
selective excitations of the triplet states shown in Fig. || 
p4j . Let us see the excitation in detail. We consider the 
ladder case shown in Fig. |I]as an example, and follow the 
description in Ref. |2^. The effective spin is introduced 
within the subspace consisting of the ground state | 00) 
and the lowest excited state | 1 — 1). This treatment 
is helpful to visualize the transition as a rotation of the 
spin. The Hamiltonians which govern the system are, 

H = Ho + Hx, (1) 
H = J ^ Sji • s i2 + gus ^ Hfo) ■ (sq + s i2 ), (2) 



Sn ■ S,i 



nn nn 



(3) 



Here, J and J\ are the intra- (rung) and inter-pair (leg) 
exchange interactions, respectively. H(^j) is the external 
field at the i-th site and "nn" means the nearest neighbor 
sites. By introducing the operators Si = Sji + and 
Tj = sn — Sja, Eqs. (|l])~(||) are rewritten as, 



ff o = \jJ2 s *+wbJ2 h ( x >)- s * ( 4 ) 

i i 

Hi = ijx^Si-Sj+Ti-Tj), (5) 



where some constant terms are dropped. 

We further rewrite these Hamiltonians in terms of the 
Pauli spin components in the following manner. 



S? = 1 - Sf = S? = 0, S* = ±«-l), 



(6) 



As a consequence of these transformations, the two states 
of | 00) and | 1 — 1) can be treated as the two spin states 
represented by the Pauli spin matrices. The total Hamil- 
tonian, Eq. ([j]) is expressed by these spin matrices as, 
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This describes the system of the spins coupled through 
anisotropic exchange interactions. The microwave ro- 
tates these spins. 

The small-angle rotations of the spins in a small re- 
gion create a packet of k ~ magnons like a soliton. 
The packet is localized on the chain due to the magnetic 
field gradient. A mismatch in the magnon excitation en- 
ergies between adjacent regions along the chain prohibits 
the packet from moving to the lower field region. On 
the other hand, the continuum excitations near the one- 
magnon excitations at k=0 is absent, so that it is dif- 
ficult for thepacket to move to the higher field region 
p7| , p9| , [30| , [bi] unless the process of the energy re- 
lease by phonon emissions is considerable. Consequently, 
the magnons are confined in the region where they are ex- 
cited, and the SN interaction is produced only between 
the qubit pair of interest. 

The population of the magnons with k = (= n(0)) is 
determined by the balance between excitation and relax- 
ation, 



dn(0) 
dt 



n(0) 
~t7' 



(8) 



where W eK is the transition probability of | 00) — >| 1 — 1) 
per unit time by the microwave irradiation, and T s is the 
magnon lifetime. At the steady state, dn(0)/dt — 0, so 
that n(0) = W ey T s . As the microwave irradiation is shut 
off, the spins start to relax to the ground state (a z = —1) 
with the relaxation time given by T s . 

We end up this section with some remarks on the na- 
ture of the excited states. One is about the lifetime of 
the excited states (T s ); it could be rather long because 
the excitation is primarily forbidden for the usual elec- 
tric dipolar transition. Although the forbiddance itself is 
usually lifted by some additional interactions such as the 
Dzyaloshinsky-Moriya interaction 24 p5| , the primarily 



forbidden transition is expected to lead to rather long 
lifetime of the excited triplet states and narrower transi- 
tion line, which is favorable for the selective excitation. 
Another is about the wave number of the triplet states; 
they should have k=0 and the transition to the staggered 
component (k — a/ir) is forbidden, because only the tran- 
sitions with the momentum transfer q = k — k =0 are 
allowed by the microwave, and the ground state ( 1 00) ) is 
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a uniform singlet state (fc=0). This fact is quite favor- 
able for the long-range inter- nuclear coupling with q ~ 0. 
Note that the magnons excited by the microwave are in 
the non-equilibrium states far from the thermal equilib- 
rium where the k = a/ir magnons are primarily excited. 



IV. INTER-NUCLEAR COUPLINGS 
MEDIATED BY MAGNONS 



We next look into the details of the longitudinal com- 
ponent of the SN interaction caused by a packet of spin 
triplet states. For simplicity, we assume the following on- 
site anisotropic hyperfme Hamiltonian in the z-th site, 

H hf = {A^sfJ! + \At{s+ir + s -/+)}, (9) 

which can be rewritten using the spin introduced in Eq. 



Hhf = \A\\ {Ifot - It) + -^.(cr+Jf + aril 



(10) 



Hence, besides the term, —^Anl*, which can be incorpo- 
rated into the Zeeman term in the nuclear Hamiltonian, 
the nuclear interaction with the spin a can be expressed 
by the anisotropic hyperfme interaction. The first term 
creates the shifts at the nuclear sites corresponding to 
H tr and the SN interactions, and the other terms give 



rise to the spin-lattice relaxation |32 



The ef- 
fect of the spin lattice relaxation is discussed in the next 
section. 

Since the transverse component of the SN interaction 
due to o ± vanishes in the field gradient, we can restrict 
ourselves to the longitudinal component. The longitudi- 
nal component of the SN interaction is given by, 



H SN = Walfll 



(11) 



where, 



JnA\\ 
N 



E 



cos{(fc — k )r,ij}. 



(12) 



Here, and e/t are, respectively, the population and the 
energy of the magnon with the wave number k, and 
is the distance between the two nuclei of interest. In the 
equilibrium states, nk is given by the Bose function for 
the given temperature. 

Actually, eq. ( pd| ) is a special case of the general for- 
mulas for the indirect spin-spin interaction, 



H md = *(r«)////, 
with the range function, 

Hnj) = J2 X(«) exp(^r. y ). 



(13) 



(14) 



Here, x{l) is the zero energy component of the general- 
ized susceptibility, and q — k — k . In the present case, 
the electronic state is not in the equilibrium state, so 
that x(<?) is different from that in the thermal equilib- 
rium. Note in particular, that the magnons with the 
wave numbers other than k ~ 0, such as k = Tt/a, are 
not excited here. 

The range to which this interaction reaches depends 
on the range function, which is determined by the form 
of x(?) as a function of q. For the pair of nuclei far from 
each other, the most important interaction comes from 
the uniform part of the susceptibility, i.e. q ~ caused 
by the scatterings of the magnons with k ~ 0. 

In the case of the transverse component of the SN inter- 
action in a 3-D system, the range function can be calcu- 
lated only from the magnon dispersions and has the form 
~ a/r cxp(— r/3a) (a is a lattice constant) |14| |. In our 
case, however, the situation is completely different from 
this case, and the interaction reaches rather long distance 
because of the following reasons. Firstly, in the longitu- 
dinal component, the number of excited magnons, as well 
as the magnon dispersions, are responsible for the range 
function (see Eq. (|l^)). Secondly, the system is 1-D so 
that a qualitative difference exists in the structures of 
x(q)- Thirdly, the system is not in the equilibrium state. 
Since x(q) is enhanced at q ~ in the present case, a 
rather long distance interaction is expected between the 
nuclei of interest. Note again that the contribution from 
the magnons with k ~ ir/a to x(q ~ 0), which exists in 
the equilibrium state, does not exist in the present case. 

Here, we make a rough estimation of the range function 
for the case of Fig. |i| using Eq. (|l2|). The magnon 
disp ersion of the triplet state in this situation is given by 
§5|, 



e(k n ) = C + J(jx - ~jf) cos(fc n ) + 



(15) 



where k n = mr/N, j\ = J\/ J and C is the part inde- 
pendent of k n . Here, the magnetic field is assumed to be 
uniform in this region. Using Eq. (15) and recalling that 
the microwave irradiation excites only the k — compo- 
nent of the magnons, i.e., n(k) = for kj^ 0, one can 
calculate the range function Wij in Eq. (|12j), 



InAn 



2 N 

E 



2n(0) 



N j ^ e(fcn) _ e(0 )- s ^) 

2 7 ^M0)/iV} ^ cos(fc„r y ) 



J{h-\jl)N ^cos(fc n )-l' 



(16) 



Assuming N = 20, r y = 10 (in units of a), A\\ = 100 
kOe//is, 7 n /(27r) = 4.3 MHz/kOe ( : H as an example), 
J = 50 K, j 2 = 0.2 and n(0)/N = 0.01, one obtains 
Wij = 15 kHz, which is the same order of magnitude 
as the nuclear dipole coupling acting between nuclei 3 
A apart from each other 0, a, and one to three orders 
of magnitude greater than the J-couplings used in the 
solution NMR-QC's EL B. 
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The estimation here is for the ideal case, because the 
population difference may be underestimated (n(k ^ 0) 
may be populated somewhat in the actual situation). 
Moreover, the realistic magnon dispersions as well as the 
conditions of the microwave excitation such as n(0)/N 
and iV may affect the value of Wij. Nevertheless, this 
estimation indicates that the longitudinal component of 
the SN interaction could be large enough to serve as a 
c-NOT gate. Note that the strength of the coupling, Wij 
can be controlled by the microwave intensity via n(0), be- 
cause it is determined by the balance between excitation 
and relaxation, and in the steady state, n(0) = W ex T s . 



V. DECOHERENCE 

Decoherence is one of the serious concerns in this 
model. On the one hand, the excited triplet states can 
create inter-qubit couplings, but on the other hand, they 
inevitably activate the scattering channels of the nu- 
clear spins and reduce the spin-lattice relaxation time 
T\ H^, |33|, Although this is unavoidable, it is still 
possible to reduce the chances of decoherence. 

Let us see the scattering process of magnons by a 
nucleus in detail. Since the direct one-magnon process 
is prohibited because of the mismatch between nuclear 
and electron excitation energies, the possible lowest or- 
der process is the two-magnon (Raman) process. This 
process can be the major cause of T\ relaxation when the 
number of excited magnons is increased. The process is 
known, however, to require highly anisotropic exchange 
interactions because of the conservation of the angular 
momentum before and after the scattering. Hence, the 
process can be reduced in the system with isotropic ex- 
change interactions |32], [33). The three-magnon process 
is the next possible process, where two incident magnons 
are scattered by a nucleus and one magnon is emitted. 
This process is active even in the case of isotropic ex- 
change interactions. 

In general, the spin-lattice relaxation time is given by, 



1 



2^A 2 ± k B T ^Im X + -(q,w„) 



E 



(17) 



here, x" 1 (q,Ll>) is the transverse component of the dy- 
namical susceptibility and u) n is the NMR frequency. 
(Note again, that the present situation is not a ther- 
mal equilibrium state, so that x(<7j w ) is different from 
those in the thermal equilibrium.) As seen in Eq. (17), 
1/Ti is proportional to A\. Recalling that An is the 
only necessary component of the hyperfine coupling for 
the longitudinal component of the SN interaction, the 
highly anisotropic hyperfine coupling (Ay /Ax 3> 1) may 
be used to reduce 1 /Ti without reducing the longitudinal 
component of the SN interaction. Such highly anisotropic 
hyperfine interaction can be realized by transfered hyper- 
fine couplings |3(| . 

Even with such effects, however, it might be inevitable 
to invoke the quantum error correction method for the 



unavoidable nuclear spin flips in the end p7[ pq . Even so, 
the chain structure is fortunate for the error correction 
process. The error correction requires three equivalent 
QC's, which are entangled like ana^ais at each qubit 
I, where an is an eigenstate (e.g. spin-up state) at the 
qubit / of the i-th QC. Such set of nuclei can be provided 
by the three consecutive chains. It is possible to align 
the qubit nuclei (see § VI) and once they are aligned, the 
inter-chain nuclear dipole couplings keep them entangled. 
(For this purpose, only the I Z I Z component is needed. 
Since I + I~ components make the NMR line broad, they 
should be eliminated by applying a small field gradient 
in this direction.) 



VI. NUCLEAR ALIGNMENTS 

So far, we have described the model for a single QC. 
Unfortunately, the sensitivity of the conventional NMR 
method is not high enough to detect signals from a single 
QC, so that one needs to integrate many equivalent QC's 
to obtain the results of computations. In this process, 
the nuclear alignment can help us to reduce the number 
of equivalent nuclei required for each qubit. It is also 
advantageous for obtaining a narrower NMR line. 

The nuclear alignment is also associated with the 
qubit initialization. Although the qubit initialization 
may be achieved by the pseudo-pure state technique J3S| ] 
and/or the algorithmic cooling |10| |d|, aligning nuclear 
spins, even partially, can assist the methods to work in 
the systems with the large number of qubits (N), be- 
cause the number of QC's which happen to be in the 
pure state in the thermal equilibrium is proportional to 
(Tiuj n /2kBT)N/2 N , which becomes smaller and smaller 
as increasing N. The nuclear alignment can relieve this 
problem. 

Here, we propose a possible configuration of a "qubit 
initializer", which partially aligns the qubits in the spin 
chains by the optical pumping technique p2[ . The 
scheme enables us to separate the materials responsible 
for the initialization and computation, so that one can 
optimize the computation part without being restricted 
by the initialization scheme. The schematic illustration 
is shown in Fig. ||. Films made of spin chains are sand- 
wiched by Si single crystals and placed in a field gra- 
dient and at a low temperature. The electrons in the 
semiconducting Si matrix are polarized by a circularly 
polarized near-infrared laser light with an energy corre- 
sponding to the semiconducting gap of Si ^] , and 
the polarizations of the electrons are transfered to the 
nuclear spin systems of th e sp in chains through the hy- 
perfine couplings |4(], ^7J, or the cross-polarization 
/coherent transfer techniques Here, the films and 

the Si matrix are not necessarily bonded chemically but 
contacted mechanically, because the direct dipolar cou- 
plings are available for the polarization transfer process 

After the nuclear alignment is completed, the laser 
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Si matrix Spin chain layer 




Circularly polarized light for optical pumping of Si 
Microwave irradiation for spin chains 



FIG. 5: Schematic illustration of a qubtt initializer for an in- 
tegrated NMR-QC and a possible experimental configuration. 
Films made of spin chains are embedded in a Si matrix and 
placed in a field gradient and at a low temperature. Small 
balls represent the qubit nuclei in the spin chains. A laser 
light for the optical pumping of the electrons in the Si matrix 
and the microwave irradiation for the spin chains are delivered 
from the bottom (Faraday configuration). 

light is shut off. The excited electrons in the Si ma- 
trix will die out quickly and the Si matrix will return 
to the silent environment. Here, the Si matrix is se- 
lected because the abundant nuclei ( 28 Si) have no spins 
(1 = 0). Since the optical pumping is effective in most of 
the semiconductors f44| . other semiconductors with no 
nuclear spins can be used as well. For example, CdTe 
enriched by the isotopes with no nuclear spins, such as 

110,112,114 Cd and 126,128,130 Te j g ^ availabl ^ 



VII. PRACTICAL ISSUES IN THE 
IMPLEMENTATION 

We have shown an ideal model for the NMR QC so far, 
but one might notice some practical issues to implement 



this model. First of all, one should find a suitable spin 
chain system with a singlet ground state, containing more 
than two stable isotopes for one element, one of which has 
/ = 1/2. Moreover, the periodical placements of the 1/2 
nuclear spins may require some ingenious techniques. 

One of the possible candidates for the QC would be 
a film made of some organic materials with spin gaps. 
They contain *H and 13 C having I = 1/2. They also 
have large unit cells, which are useful to reserve large 
spacial distances between qubits. Moreover, one could 
utilize a well-established chemical technique of the selec- 
tive isotope replacement of 1 H by 2 D (7 = 1) Q or 12 C 
(I = 0) by 13 C (I = 1/2). The periodic placements of 
the qubits may be achieved with the epitaxial growth pc| ] 
and/or the Langmuir-Blodgett methods as follows. One 
prepares two sets of molecules with the same chemical 
formulas but different isotopes by the selective isotope 
replacement technique, for example, all deuterized sam- 
ples and those having 1 H's at one of the hydrogen sites. 
These two sets of molecules could be layered epitaxially 
using the MBE or the Langmuir-Blodgett method so as 
for -"^H's or 13 C's to be placed periodically. 

In order to avoid unnecessary nuclear couplings, the 
nuclei other than qubits are preferable to have no spins 
(1 = 0). In this sense, it is fortunate that major abundant 
isotopes in organic materials such as 12 C and 16 have no 
spins. The Haldane systems such as NENP and NINO 
are also fortunate because they contain 58 Ni, an abun- 
dant isotope with no spins. In addition, a high power 
decoupler and/or a decoupling sequence between unlikc- 
spins using a train of 7r-pulses can counteract the effects 
of the nuclei with 1^0 such as 2 D. 
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